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The  technologies  and  the  processes  for  the  use  of  biomass  as  an  energy  source  are  not  always  environ¬ 
mental  friendly.  It  is  worth  to  develop  approaches  aimed  at  a  more  sustainable  exploitation  of  biomass, 
avoiding  whenever  possible  direct  combustion  and  rather  pursuing  fuel  upgrade  paths,  also  considering 
direct  conversion  to  electricity  through  fuel  cells.  In  this  context,  it  is  of  particular  interest  the  develop¬ 
ment  of  the  biomass  gasification  technology  for  synthesis  gas  (i.e.,  syngas)  production,  and  the  utilization 
of  the  obtained  gas  in  fuel  cells  systems,  in  order  to  generate  energy  from  renewable  resources.  Among 
the  different  kind  of  fuel  cells,  SOFCs  (solid  oxide  fuel  cells),  which  can  be  fed  with  different  type  of  fuels, 
seem  to  be  also  suitable  for  this  type  of  gaseous  fuel.  In  this  work,  the  syngas  composition  produced  by 
means  of  a  continuous  biomass  steam  gasifier  (fixed  bed)  has  been  characterized.  The  hydrogen  concen¬ 
tration  in  the  syngas  is  around  60%.  The  system  is  equipped  with  a  catalytic  filter  for  syngas  purification 
and  some  preliminary  tests  coupling  the  system  with  a  SOFCs  stack  are  shown.  The  data  on  the  syn¬ 
gas  composition  and  temperature  profile  measured  during  the  experimental  activity  have  been  used  to 
calibrate  a  2-dimensional  thermodynamic  equilibrium  model. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  gas  produced  through  gasification  process  (i.e.,  syngas)  can 
be  used  in  several  applications  for  separate  or  combined  heat  and 
power  generation,  and  for  the  production  of  liquid  fuels  and  chem¬ 
icals.  Several  goals  have  been  achieved  in  the  last  20  years  in  this 
research  field,  but  the  process  efficiency  and  the  tar  removal  are 
still  a  problem  for  most  of  the  applications. 

The  syngas  produced  can  be  used  for  power  generation  in  a 
CHP  cycle  (combined  cycle  for  heat  and  power  generation).  In 
these  systems  the  gas  yield  can  be  exploited  in  alternative  inter¬ 
nal  (IC)  combustion  engines.  Modified  gas  engines  can  be  fed  with 
gas  of  different  quality,  even  those  having  calorific  values  around 
5MJNm-3.  The  output  energy  can  be  roughly  divided  up  as  one- 
third  electric  energy  and  two-third  thermal  energy. 

IGCC  (integrated  gasification  combined  cycles)  are  preferred  - 
instead  of  CHP  systems  based  on  alternative  engines  -  for  the  elec¬ 
tricity  generation  on  larger  scale.  In  these  systems  the  synthesis  gas 
is  burned  in  the  combustion  chamber  of  a  gas  turbine.  However,  gas 
turbines  requires  a  pressurized  gas  supply  line,  thus  the  syngas  has 
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to  be  compressed  or  directly  produced  in  a  pressurized  gasification 
unit  at  a  pressure  of  5-20  bar.  The  major  drawback  for  pressur¬ 
ized  gasification  is  the  cost  related  to  the  energy  consumption  to 
pressurize  the  gasifying  agent  and  the  inlet  biomass. 

Considering  the  existing  experimental  pilot  plants,  fuel  cells  can 
potentially  achieve  higher  electrical  efficiencies,  compared  with 
conventional  internal  combustion  engines.  A  fuel  cell  installed 
downstream  of  a  gasifier  -  after  an  enhanced  cleanup  stage  -  repre¬ 
sents  a  coupled  system  (similar  in  size  to  an  equivalent  gasification 
-  micro  gas  turbine  generation  system)  characterized  by  high  elec¬ 
tric  and  overall  efficiencies. 

As  confirmed  by  previous  research  studies  available  in  the  lit¬ 
erature,  the  use  of  synthesis  gas  from  biomass  gasification  in  fuel 
cells  is  still  in  an  early  stage  of  application.  The  performance  of  inte¬ 
grated  biomass  gasifiers  and  SOFC  systems  has  been  investigated 
from  a  theoretical  point  of  view  by  several  authors  [1-9]. 

Panopoulos  et  al.  [2,3]  studied  the  feasibility  of  high  efficiency 
SOFC-CHP  systems  of  sizes  up  to  1  MWe  with  a  novel  allother- 
mal  steam  gasification  reactor.  The  reactor  is  heated  by  burning 
depleted  SOFC  off  gas,  un-reacted  residual  char  and  additional 
biomass  if  necessary.  The  STBR  used  is  0.6  to  limit  the  heat  required 
for  the  gasification  system.  The  size  of  the  system  is  based  on 
100  m2  of  SOFC  active  surface.  The  model  results  show  a  system 
electrical  efficiency  of  36%  and  a  thermal  efficiency  of  14%. 
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Nomenclature 

A  gasifier  cross  section  area  (m2) 

a  thermal  diffusivity  (m2  s-1 ) 

cp  specific  heat  (kj  kg-1 I<-1 ) 

H  specific  enthalpy  (kj  kg-1 ) 

AH  specific  enthalpy  variation  (kj  kg-1 ) 
m  mass  flow  rate  (kg  s-1 ) 

n  number  of  moles 

n  unit  vector  normal  to  the  relevant  surface 

N  Input  vector  in  the  thermodynamic  model 

p  pressure  (Pa) 

r  gasifier  radial  coordinate 

Q  heat  source  ( W  m-3 ) 

Q.  thermal  power  (W) 

T  temperature  (I<) 

u  velocity  vector  (ms-1) 

v  velocity  component  along  the  radial  coordinate,  r 

(ms-1) 

w  velocity  component  along  the  vertical  coordinate,  z 

(ms-1) 

W  average  velocity  on  the  gasifier  cross  section  (m  s-1 ) 
z  gasifier  vertical  coordinate 

V  biomass  volume  (m3) 

Greek 

r]  efficiency 

rjc  carbon  conversion  efficiency 

0  porosity 

X  thermal  conductivity  (W  m-1  K-1 ) 

p  density  (kg  m-3) 

y  adimensional  parameter 

§  biomass  fraction  converted  into  gas 

r  time  constant  (s) 

Subscripts 

0  values  at  the  initial  step  (at  the  gasifier  inlet) 

A  air 

B  biomass 

C  carbon 

E  heating  condition  on  the  external  surface  of  the  gasi¬ 

fier  shell 
G  syngas 

H  hydrogen 

O  oxygen 

i  values  at  the  ith  step  (in  the  gasifier) 

IN  values  at  the  gasifier  inlet 

OUT  values  at  the  gasifier  outlet 

Q  flow  (NmL  min-1) 

R  reaction 

S  steam  condition  on  the  internal  coaxial  tube  of  the 

gasifier 

VAP  vaporization 

W  water 

Acronyms 

CFD  computational  fluid  dynamic 

LHV  lower  heating  value 

PFD  process  flow  diagram 

RT  room  temperature 

SC  steam  to  carbon:  ratio  between  the  moles  of  FI20  fed 

and  the  moles  of  carbon  in  the  feedstock 
STBR  steam  to  biomass  ratio 


Cordiner  et  al.  [4]  simulated  a  14kW  energy  generation  system 
consisting  of  a  SOFC  stack  fed  by  a  biomass  gasifier.  The  gasifier  is 
integrated  with  a  heat  recovery  system  on  the  SOFC  off-gas  com¬ 
bustion,  in  order  to  enhance  the  energy  performance.  The  gasifier 
has  been  simulated  by  means  of  a  zero-dimensional  equilibrium 
model,  while  a  3D  CFD  model  has  been  implemented  for  the  SOFC. 
The  overall  system  efficiency  has  been  assessed  at  45%. 

Omosun  et  al.  [5]  have  developed  a  steady  state  model  to  investi¬ 
gate  the  performance  of  an  integrated  SOFC  stack  -  biomass  gasifier 
system.  The  model  has  been  used  to  study  the  system  efficiency  and 
its  costs  considering  two  different  options:  “cold”  gas  cleaning  at 
low  temperature  and  “hot”  gas  cleaning  at  high  temperature.  The 
simulations  show  an  overall  system  efficiency  of  60%  for  the  hot 
gas  cleanup  and  34%  for  the  cold  one. 

Liu  et  al.  [6]  have  assessed  the  theoretical  performance  of 
an  integrated  gasifier-SOFC  pilot  system  which  was  under  con¬ 
struction.  It  consists  of  a  fixed  bed  gasifier  and  a  5kW  SOFC 
CHP  system.  Two  solutions  for  gas  cleaning  have  been  con¬ 
sidered:  a  combined  low  and  high  temperature  gas  cleaning 
system  and  a  high  temperature  gas  cleaning  system.  The  results 
of  the  analysis  show  that  the  electrical  efficiency  of  the  CLIP 
system  considering  the  different  gas  cleaning  units,  are  approxi¬ 
mately  the  same,  although  the  high  temperature  solution  shows 
higher  thermal  efficiency  both  for  the  energy  and  exergy  bal¬ 
ance. 

Nagel  et  al.  [7-9]  have  researched  on  the  technical  and  on 
the  cost  analysis  of  a  biomass  integrated  gasification  fuel  cell  (B- 
IGFC)  system.  The  paper  [9]  shows  some  preliminary  experimental 
results  collected  from  a  series  of  tests  performed  in  the  B-IGFC 
demonstration  unit  consisting  of  gasifier,  gas  cleaning  session  (with 
catalyst)  and  fuel  cell  stack.  The  gasifier  and  the  catalytic  unit  could 
operate  with  a  satisfying  degree  of  reliability;  the  SOFC  unit  fed 
with  syngas  generates  approximately  40%  less  current  compared 
to  methane  operation.  Ash  deposition  was  the  major  obstacle  for  a 
smooth  SOFC  system  operation. 

The  purpose  of  the  present  work  is  the  evaluation  of  the  feasibil¬ 
ity  of  a  B-IGFC  for  residential  applications  (i.e.,  micro-cogeneration 
system).  The  main  interest  is  the  analysis  of  the  steam  gasifica¬ 
tion  process  as  a  possible  pathway  to  generate  high  quality  syngas 
suitable  for  solid  oxide  fuel  cells. 


2.  Modeling 

2.1.  Process  conditions 

The  aim  of  the  present  work  is  a  feasibility  assessment  of  a 
B-IGFC  micro-cogeneration  system  suitable  for  residential  applica¬ 
tions.  The  electric  energy  load  required  by  1  -2  apartments  has  been 
considered  for  choosing  the  output  of  the  system  (3-5  kWel ).  There¬ 
fore,  the  gasification  unit  has  been  designed  considering  a  SOFCs 
stack  able  to  provide  5  kWel.  The  gasifier  input  power  is  between 
11  and  12kWth,  estimated  considering  a  SOFCs  stack  efficiency  of 
45%  and  also  taking  into  account  the  thermal  losses  through  the 
gasification  system.  The  test  feedstock  (wood  pellets)  has  a  lower 
heating  value  of  1 8-1 9  MJ  kg-1 . 

The  characterization  of  the  gasification  process  has  been  carried 
out  considering  the  fuel  gas  required  by  the  solid  oxide  fuel  cells. 
Fuel  gas  for  SOFCs  has  to  be  particle  free,  in  order  to  avoid  clogging 
of  the  gas  channels,  sulphur  compounds  free  to  avoid  the  poisoning 
of  the  cells  and  must  have  low  tar  content.  In  fact,  the  cells  can 
tolerate  a  certain  amount  of  tars  because  they  can  directly  reform 
hydrocarbons  which  is  also  a  way  for  the  chemical  cooling  of  the 
SOFC  itself.  This  type  of  fuel  cells,  which  work  at  high  temperature, 
can  accept  different  fuels  in  input,  even  if  the  standard  fuels  are 
hydrogen  and  methane. 
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Table  1 

Syngas  molar  composition  (dry  basis)  and  char  production,  LHV  and  HHV  and  global  process  efficiencies. 


(%  molar 
fraction) 

Partial  oxidation 

ER  =  0.20,  7=740  °C 

Steam  gasification 

SC  =  0.50,  7=720°C 

Steam  gasification 

SC  =  3.00,  7  =  750  °C 

h2 

29.35 

51.93 

62.80 

CO 

33.90 

37.16 

13.87 

co2 

5.23 

8.36 

23.26 

ch4 

0.62 

2.50 

0.03 

n2 

30.90 

0.05 

0.04 

02 

0.00 

0.00 

0.00 

Char 

0.00 

0.00 

0.00 

LHV  (MJ  Nirr3) 

7.57 

11.05 

8.43 

HHV  (MJNirr3) 

8.16 

12.16 

9.65 

Global  process  efficiency  (%) 

86.4 

83.2 

41.0 

Syngas  generated  through  biomass  gasification  shows  higher 
hydrogen  concentrations  when  pure  steam  or  a  mixture  of 
air-steam  is  used  as  gasifying  agent.  Taking  into  account  the  need 
to  maximize  the  hydrogen  concentration  in  the  syngas,  the  pure 
steam  gasification  process  has  been  considered  as  the  most  suitable 
process  for  the  purpose  of  this  work. 

The  test  gasifier  is  a  downdraft  fixed  bed  reactor  with  a  cylindri¬ 
cal  shape  and  it  is  externally  heated  by  means  of  4  electric  ovens. 
The  steam  is  provided  by  an  internal  coaxial  tube  (see  Section  3.1 ). 
This  particular  configuration  has  been  chosen  to  maximize  the  heat 
exchanges  between  the  steam  -  fed  at  high  temperature  (up  to 
600  °C)  -  and  the  biomass  bed  in  the  reactor.  The  reactor  geometry 
has  been  designed  to  ensure  the  nominal  load  capacity  (2.5  kg  h- 1 ) 
required  to  provide  the  11-13  kWth  input  load. 

2.2.  Thermodynamic  approach 

A  thermodynamic  analysis  has  been  carried  out  by  means  of  a 
2-phase  (gas-solid)  thermodynamic  model  applied  to  the  biomass 
conversion  stage,  based  on  the  Cantera  libraries  [10].  The  imple¬ 
mented  solution  algorithm  is  a  version  of  the  Villars-Cruise-Smith, 
based  on  the  Gibbs  energy  minimization  method.  The  model  takes 
into  account  60  different  species  for  the  gas  phase,  and  1  for  the 
mass  phase  (i.e.,  graphite)  [11].  The  NASA  [12]  and  Gri-Mech  [13] 
databases  have  been  used  for  the  gas,  water  steam  and  char  ther¬ 
modynamic  properties. 

From  the  equilibrium  analysis,  and  for  the  considered  feed¬ 
stock  (i.e.,  spruce  pellets),  the  optimum  process  efficiency  has  been 
obtained  by  means  of  partial  oxidation.  However,  from  a  practical 
point  of  view,  the  steam  gasification  gives  the  highest  syngas  qual¬ 
ity,  since  it  is  characterized  by  a  higher  concentration  of  hydrogen 
and  it  is  not  diluted  by  the  atmospheric  nitrogen  as  in  the  partial 
oxidation  case;  moreover,  there  is  a  lower  production  of  char,  due 
to  a  higher  carbon  conversion  efficiency.  On  the  other  hand  partial 
oxidation  is  a  simpler  conversion  process  as  it  has  usually  less  prac¬ 
tical  heat  transfer  limitations  for  high  temperature  heat  exchanges 
than  steam  gasification.  The  global  process  efficiency  of  both  partial 
oxidation  and  steam  gasification  have  been  computed  taking  into 


account  the  feedstock  and  the  syngas  heating  values,  the  reactions 
enthalpies  and  the  enthalpies  of  the  gasifying  agents  (1 ). 

_  LHVg  ■  riiG  —  AHr  —  AH  a  ■  riiA  —  ( AHw  —  A  Hvap  —  A  Hs )  ■  rilw  ,  -  >. 

11  ~  LHVb  ■  ifiB 

Tables  1  and  2  contain  the  results  obtained  for  partial  oxida¬ 
tion  process  with  ER  =  0.20,  740  °C  conversion  temperature  and  for 
steam  gasification  with  SC  =  0.5  and  720  °C  conversion  temperature. 
These  values  of  the  process  parameters  determine  the  operation 
conditions  that  lead  to  the  optimum  global  efficiencies  of  the  con¬ 
sidered  gasification  process.  Furthermore,  it  has  also  considered  the 
steam  gasification  process  with  SC  =  3  and  750  °C,  that  approaches 
to  more  realistic  plant  operation  conditions. 

2.3.  FEM  approach 

The  thermodynamic  analysis  of  the  designed  steam  gasifica¬ 
tion  reactor  has  been  extended  on  a  2-dimensional  domain  of  the 
designed  gasifier,  coupling  the  heat  and  the  momentum  equations 
and  solving  the  PDE’s  problem  by  means  of  a  finite  elements  for¬ 
mulation. 

The  proposed  model  is  based  on  a  system  of  two  PDE’s,  the  heat 
conduction  Eq.  (2)  and  the  Euler  Eq.  (3),  defined  in  the  stationary 
regime  and  in  two  dimensions.  The  flow  field  of  the  gasifier  bed 
has  been  characterized  by  means  of  the  Euler  equation,  since  the 
viscous  and  volume  forces  have  been  considered  negligible  for  the 
present  analysis.  For  the  purposes  of  the  present  analysis  the  bed 
material  can  be  considered  as  a  granular  fluid  and  the  resulting 
average  value  of  its  velocity  is  very  low.  The  coupling  of  Eqs.  (2) 
and  (3)  gives  the  temperature  T  and  the  velocity  u  (4)  distributions 
as  results. 

(XV2T)  = -Q(T)  (2) 

p(u  •  V)u  =  -Vp  (3) 

u  =  (uy,  uz)  =  (i/,w)  (4) 

The  PDE’s  discretization  procedure  has  been  performed  with  the 
finite  element  analysis  by  means  of  the  Comsol  Multiphysics®  soft¬ 
ware  [14].  The  continuous  computational  domain  (i.e.,  the  biomass 


Table  2 

Enthalpy  variation  along  the  conversion  process  for  several  values  of  the  gasifying  agent  inlet  temperature. 


Enthalpy  variation  (kWh  kg-1 ) 

Partial  oxidation 

ER  =  0.20,  7=740°C 

Steam  gasification 

SC  =  0.50,  7=  720°C 

Steam  gasification 

SC  =  3.00,  7=750°C 

^gasifying  =  20  °C 

0.58 

^gasifying  =  1 00  °  0 

0.55 

1.42 

1.98 

^gasifying  =  200  °C 

0.52 

1.41 

1.94 

^gasifying  =  300  °C 

0.49 

1.39 

1.83 

^gasifying  =  400  °  C 

0.46 

1.37 

1.71 

^gasifying  =  500  °C 

0.43 

1.35 

1.59 

^gasifying  =  000  °  C 

0.39 

1.33 

1.47 
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Fig.  1.  2D-axisymmetric  reactor  geometry. 


bed  and  the  reactor  walls)  has  been  divided  up  into  triangular  ele¬ 
ments,  so  that  the  unknown  temperature  field  could  be  expressed 
by  means  of  shape  functions.  These  polynomial  functions  are  com¬ 
pletely  defined  inside  each  element  by  the  temperature  value  in 
the  nodal  points.  Thanks  to  the  cylindrical  shape  of  the  reactor,  the 
whole  system  has  been  discretized  by  a  2D-axisymmetric  geom¬ 
etry  (Fig.  1).  The  symmetry  is  characterized  by  a  cross  sectional 
plane,  perpendicular  to  the  upper  and  lower  faces  of  the  cylindrical 
reactor  and  delimited  by  the  symmetry  axis. 

As  boundary  conditions  for  the  heat  conduction  equation,  tem¬ 
perature  values  have  been  assigned  to  the  reactor  external  wall  (TE) 
and  to  the  internal  coaxial  tube  side  (T5)  where  the  steam  flows. 
Two  adiabatic  conditions  have  been  assumed  for  the  inlet  (upper 
side)  and  outlet  sections  (lower  side)  of  the  reactor  (5). 

n  •  (A.VT)  =  0  (5) 


The  steam  temperature  (Ts)  has  been  estimated  by  means  of  an 
iterative  algorithm.  As  a  first  guess,  the  steam  temperature  in  the 
internal  coaxial  tube  is  set  at  the  inlet  value  (temperature  at  the 
super-heater  outlet).  By  means  of  several  iterations,  it  is  possible 
to  compute  the  heat  fluxes  entering  (or  exiting)  the  steam  tube, 
which  are  in  equilibrium  with  the  vertical  temperature  profile, 
which  becomes  the  actual  boundary  condition. 

The  boundary  conditions  assumed  for  the  Euler  equation  are 
represented  by  slip  conditions  on  the  reactor  walls,  assigned  pres¬ 
sure  (atmospheric  value)  at  the  reactor  outlet  and  assigned  velocity 
(i.e.,  average  value  on  the  cross  section  A)  (6)  at  the  reactor  inlet 


W  = 


inB 
Pb  -  A 


(6) 


The  model  is  structured  by  means  of  a  three  steps  sequence.  In 
the  first  step  the  source  term  of  the  heat  equation  is  computed  using 
a  thermodynamic  equilibrium  routine  as  a  function  of  the  process 
parameters  (temperature,  pressure,  steam  to  carbon  ratio).  Then 
the  temperature  and  the  velocity  (plug  flow  regime)  fields  are  eval¬ 
uated  with  a  finite  element  formulation  in  the  second  step;  in  this 
phase  a  2D-axisymmetric  representation  of  the  reactor  has  been 
considered.  The  final  computational  step  is  the  calculation  of  the 
reaction  products  profiles  in  the  reactor  in  equilibrium  with  the 
computed  temperature  field.  The  characterization  of  the  syngas  is 
carried  out  estimating  the  local  value  of  the  concentration  of  the 
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Fig.  2.  Enthalpy  variation  along  the  process. 


main  gaseous  compounds  (H2,  CO,  CH4,  C02,  H20)  and  comput¬ 
ing  the  average  heating  value.  Moreover,  the  solid  product  residual 
fraction  (char)  -  in  equilibrium  with  the  gaseous  one  -  is  also  cal¬ 
culated. 

The  first  computational  step  gives  the  heat  equation  source  term 
as  a  function  of  the  temperature,  considering  its  dependence  on 
the  enthalpy  variation  along  the  conversion  process,  AH.  The  value 
of  the  enthalpy  variation  has  been  estimated  using  the  aforemen¬ 
tioned  thermodynamic  approach.  The  resulting  curves  (Fig.  2)  - 
representing  the  enthalpy  variations  as  a  function  of  the  tempera¬ 
ture  -  have  been  parameterized  for  different  values  of  SC  ratio  (i.e., 
the  steam  to  carbon  ratio)  and  of  Ts  (i.e.,  the  steam  temperature). 
The  curves  have  been  fitted  with  polynomial  functions. 

The  representation  of  the  source  term  (7)  of  the  heat  equation  is 
based  on  the  estimated  polynomials  and  on  the  biomass  flow  rate 
that  is  subjected  to  the  thermal  treatment. 

Q(T)  =  —A H(T)  •  ^  =  — [Hout(T)  -  HIN(T)1  ■  ^  (7) 

The  simulation  has  been  carried  out  dividing  the  reactor  in  sev¬ 
eral  elements  (i=  1,2,...,  N)  along  the  vertical  coordinate  (z)  and 
applying  the  FEM  procedure  to  each  subdomain  step  by  step  (Fig.  3). 
For  the  present  simulations,  the  total  number  of  subdomains  has 
been  assumed  as  N  =  1 0  (i.e.,  A z  =  48  mm).  The  net  available  thermal 
power  is  calculated  as  the  difference  between  the  inlet  (Qin)  and 
outlet  (  Qout )  thermal  power.  The  design  value  of  the  inlet  thermal 
power  has  been  assumed  at  1 1  -1 2  kW. 

A  correction  factor  (§)  has  been  computed  and  applied  to  the 
heat  source  term  (8),  for  each  reactor  subdomain  next  to  the  pre¬ 
vious  one  (i.e.,  generic  ith  element).  This  factor  (9)  represents  the 
biomass  fraction  converted  along  the  gasification  process  down  to 
the  ith  step. 


Q,+1(T)  =  -A 

(8) 

e  (Qin  -  QoutV< 

^k=H\/R)f*  AHk(T,r)dr 

(9) 
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Fig.  3.  Mass  and  energy  balances. 


The  biomass  fraction  that  is  gasified  depends  on  the  ratio 
between  the  available  thermal  power  (supplied  to  the  system)  and 
the  energy  required  to  sustain  the  gasification  reactions  (i.e.,  the 
enthalpy  variation  along  the  process).  The  former  is  computed  as 
difference  between  the  inlet  and  the  outlet  thermal  flux,  while  the 
latter  is  computed  as  the  average  value  on  the  reactor  cross  section. 

Furthermore,  considering  that  the  initial  mass  of  feedstock 
entering  the  system  is  continuously  converted  into  syngas  and 
char,  the  bulk  physical  properties  of  the  biomass  bed  vary 
along  the  z-axis  (10)-(13).  The  pellets  porosity,  whose  aver¬ 
age  initial  value  (0/N  =  0.49)  is  deduced  from  the  particle  density 
(11 50-1400  kg  m-3)  and  the  bulk  density  (600-650  kg  nrr3)  of  the 
material  as  stated  by  the  producer,  have  been  taken  as  reference 
parameters. 


0i  =  (1  +§)-0/N 

(10) 

A+l  =  4>i  •  Pc  +  0  -  0f)  •  PB 

(ii) 

cPi+ 1  =  4>i  ■  CPG  +  (1  +  4>i)  •  CpB 

(12) 

^i+ 1  =  4>i  •  +  (1  —  0i)  • 

(13) 

The  last  step  of  the  analysis  is  represented  by  the  syngas  charac¬ 
terization.  For  this  purpose  an  estimation  of  the  radial  distribution 
of  the  main  chemical  compounds  present  in  the  gas  phase  and  of  the 
char  production  has  been  carried  out.  The  results  of  the  thermody¬ 
namic  equilibrium  computation  are  given  as  H2,  CO,  C02,  CFI4,  FI20 
concentrations  and  as  char  fraction  produced,  so  that  it  is  possible  to 
compute  the  process  conversion  efficiency.  For  each  computational 
step  (i.e.,  subdomain),  the  temperature  radial  distribution  obtained 
from  the  PDEs  problem  have  been  used  to  compute  a  local  estimate 
(r,  z)  of  the  equilibrium  products  concentration  in  the  reactor  bed 
(Fig.  4). 

3.  Experimental  procedures 

The  experimental  apparatus  designed  and  built  is  a  CHP  system 
based  on  a  5kWel  SOFC  stack  (SOFCpower)  fueled  with  biomass 
syngas. 

3. 1.  Gasification  CHP  unit 

The  steam  gasification  process  is  performed  at  atmospheric 
pressure  in  a  llkWth  (inlet  power)  co-current  fixed  bed  gasi¬ 
fier,  consisting  of  two  coaxial  cylindrical  Inconel®  vessels  ( 1 1 00  °C, 


- H2 

— H20 

— 1 —  CO 
— n— C02 
— ' •— CH4 


Fig.  4.  Temperature  distribution  in  the  gasifier  and  local  estimate  of  the  equilibrium 
products  concentration  (z  =  432  mm)  in  the  reactor  (SC  =  2,  7^  =  1 073,  Ts  =  873). 


2  barg  design),  having  external  diameter  of  102  mm,  height  of 
480  mm,  and  net  available  volume  of  0.0037  m3.  An  internal  coax¬ 
ial  tube  of  smaller  diameter  provides  the  gasifying  agent  (i.e., 
superheated  steam).  The  gasifying  unit  is  equipped  with  a  steam 
generator  (evaporator  and  super-heater)  capable  to  supply  steam 
up  to  600  °C  at  2  barg. 

The  gasifier  is  externally  heated  by  means  of  four  ovens  (max¬ 
imum  power  1  kW  each).  Temperature  and  pressure  sensors  are 
placed  inside  the  reactor  in  order  to  control  the  process.  A  power 
meter  measures  the  energy  consumption  of  the  ovens  and  of  the 
steam  generator  during  the  heating  phase  and  the  gasification  stage 
(stationary  temperature  stage).  The  biomass  is  stored  in  a  closed 
hopper  and  fed  in  the  reactor  from  the  top,  by  means  of  an  hori¬ 
zontal  Archimedean  screw.  The  feeding  rate  depends  on  the  rpm 
(rounds  per  minute)  of  the  screw  which  can  be  manually  set  by  an 
inverter.  The  minimum  feeding  rate  is  10  g  min-1. 

A  second  Archimedean  screw  (of  the  same  size)  is  placed  at  the 
bottom  of  the  reactor  to  discharge  the  char  produced  during  the 
gasification  process.  The  starting  and  switching  off  of  the  screws  is 
controlled  via  computer.  The  gasifier  can  work  both  in  continuous 
and  semi-continuous  mode. 

Wood  pellets  (spruce)  have  been  used  as  feedstock  for  the 
experimental  campaign  (elemental  analysis  and  heating  value  are 
reported  in  Table  3).  The  plant  has  been  tested  with  pellets  charge 
of  1,  1.5  and  2kgh-1  (up  to  80%  of  the  nominal  load  which  is 
2.5  kg  h-1). 
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Table  3 

Feedstock  properties. 

Ultimate  analysis  (%  mass  a.r.) 

C:  47.0  H:  6.1  O:  46.8  N:  0.01  S:  <0.01 

Moisture  content  (%  mass  a.r.)  7.0 

Ash  (%  mass  a.r.)  0.2 

Lower  heating  value  (MJ  kg-1)  18.5 


The  syngas  is  cooled  and  cleaned  before  performing  the  analysis 
by  means  of  a  portable  Hp  Agilent  3000  Micro  GC  able  to  measure 
the  concentration  of  C02, 02,  N2,  CH4,  C2H2,  C2H4,  C2H6,  H2  and  CO. 
The  same  GC  has  been  used  to  detect  the  hydrogen  sulphide  con¬ 
centration  changing  the  gas  carrier  from  Argon  to  Helium  (Helium 
allows  the  measurements  of  CH4,  CO,  C02,  N2  and  H2S). 

The  temperature  field  inside  the  reactor  is  characterized  by 
means  of  three  K-thermocouples  placed  at  different  height  and 
radial  positions.  The  PFD  of  the  system  is  reported  in  Fig.  5. 

3.2.  Syngas  characterization 

The  aim  of  the  experimental  activity  is  to  produce  a  syngas  suit¬ 
able  for  solid  oxide  fuel  cells,  that  implies  high  hydrogen  content 
and  low  tar  concentration.  The  experimental  conditions  to  maxi¬ 
mize  the  hydrogen  concentration  (temperature  and  SC)  have  been 
assessed  by  means  of  the  equilibrium  model  described  in  Section 
2. 

Fifteen  tests  have  been  run  during  the  experimental  activity,  as 
reported  in  Table  4.  The  gasification  temperature  investigated  has 
been  changed  between  700  and  800  °C.  The  experimental  activity 
has  been  divided  in  two  phases:  in  the  first  one  (tests  from  1  to  7th) 
the  influence  of  SC  (changed  between  2  and  3),  reaction  tempera¬ 
ture  and  steam  inlet  temperature  (from  200  °C  to  600  °C)  has  been 
investigated.  The  reactor  operated  in  a  semi-continuous  mode:  the 
biomass  has  been  fed,  and  the  char  discharged  once  every  hour. 


In  the  second  phase  (tests  from  9  to  15th)  the  attention  has 
been  focused  on  the  H2S  measurement  with  and  without  the  pres¬ 
ence  of  the  catalyst.  The  gasification  parameters,  SC  and  steam 
temperature,  have  been  kept  constant  at  2.5  and  600  °C  respec¬ 
tively.  The  ovens  temperature  has  been  kept  at  800  °C.  In  this 
way  the  experimental  average  temperature  measured  near  the 
syngas  outlet  is  about  750  °C.  This  value  represents  the  average 
conversion  temperature  and  it  is  in  agreement  with  the  conditions 
given  by  the  equilibrium  simulations  to  maximize  the  hydrogen 
concentration  and  the  process  efficiency  (cf.  Section  2.2).  In  this 
phase  the  gasifier  operated  in  a  continuous  mode:  the  biomass 
and  char  have  been  continuously  added  and  discharged,  respec¬ 
tively.  The  feeding  rate  has  been  increased  to  2  kg  h- 1 .  Test  no. 
8  has  been  also  run  in  order  to  measure  the  energy  consump¬ 
tion  of  the  empty  system.  Finally,  two  preliminary  tests  coupling 
the  system  with  the  SOFCs  stack  have  been  run  (cf.  Section 
3.4). 

By  means  of  the  known  amount  of  nitrogen  fed  to  the  sys¬ 
tem  through  the  feeding  hopper,  the  syngas  production  has  been 
estimated  to  be  around  0.6-0.7  Nm3  kg_1pellets.  The  char  produced 
during  the  steam  gasification  tests  is  approximately  18%  of  the 
initial  biomass  weight. 

In  Table  5  the  average  gas  composition  (N2  free)  has  been 
reported  for  the  12  tests  of  which  the  gas  composition  has  been 
analyzed.  No  remarkable  differences  can  be  observed  in  the  first 
series  of  tests,  in  spite  of  the  different  SC  values.  The  hydrogen  con¬ 
centration  is  higher  than  60%,  the  methane  production  is  very  low 
and  the  carbon  monoxide  and  carbon  dioxide  concentrations  have 
been  assessed  at  5%  and  30%,  respectively.  A  change  in  the  average 
gas  composition  has  been  observed  in  the  second  group  of  tests 
(9-1 2th)  due  to  the  different  operation  mode  of  the  gasifier;  indeed, 
because  of  continuous  operation,  the  gas  residence  time  inside  the 
reactor  is  reduced  and  this  implies  lower  H2  and  a  higher  CH4  and 
CO  concentrations  (approx.  53%,  7%  and  12%,  respectively).  More¬ 
over,  the  gas  heating  value  increases  in  the  second  group  of  tests 
due  to  higher  concentration  of  fuel  gases. 
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Fig.  5.  PFD  of  the  system. 
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Table  4 

Tests  performed  during  the  experimental  activity. 


Test  no. 

T  reactor 

T  steam 

SC 

Feeding 

Catalyst 

Gas  analysis 

Fuel  cell 

1 

800  °C 

200  °C 

2 

1  kg  h- 1 

- 

H2,  CO,  co2,  ch4,  c2h2,  c2h4,  c2h6 

- 

2 

800  °C 

600  °C 

2 

1  kgh-1 

- 

h2,  CO,  co2,  ch4,  c2h2,  c2h4,  c2h6 

- 

3 

800  °C 

400  °C 

3 

1.5  kg  h_1 

- 

h2,  CO,  co2,  ch4,  c2h2,  c2h4,  c2h6 

- 

4 

800  °C 

400  °C 

2 

1/1.5  kgh-1 

- 

h2,  CO,  co2,  ch4,  c2h2,  c2h4,  c2h6 

- 

56 

700  °C 

400  °C 

2 

1.5  kg  h_1 

- 

h2,  CO,  co2,  ch4,  c2h2,  c2h4,  c2h6 

- 

6 

700  °C 

400  °C 

3 

1/1.5  kgh-1 

- 

h2,  CO,  co2,  ch4,  c2h2,  c2h4,  c2h6 

- 

7 

700  °C 

400  °C 

3 

1/1.5  kgh-1 

- 

h2,  CO,  co2,  ch4,  c2h2,  c2h4,  c2h6 

- 

8 

800  °C 

600  °C 

2.5 

- 

- 

- 

- 

9 

800  °C 

600  °C 

2.5 

2  kg  h_1 

- 

h2,  CO,  co2,  ch4,  c2h2,  c2h4,  c2h6 

- 

10 

800  °C 

600  °C 

2.5 

2  kgh-1 

- 

h2,  CO,  co2,  ch4,  c2h2,  c2h4,  c2h6 

- 

11 

800  °C 

600  °C 

2.5 

2  kgh-1 

- 

h2s,  CO,  co2,ch4,c2h2,c2h6 

- 

12 

800  °C 

600  °C 

2.5 

2  kgh-1 

Yes 

h2s,  CO,  co2,ch4,c2h2,c2h6 

- 

13 

800  °C 

600  °C 

2.5 

2  kgh-1 

Yes 

h2s,  CO,  co2,  ch4,  c2h2,  c2h6 

- 

14 

800  °C 

600  °C 

2.5 

2  kgh-1 

Yes 

- 

Yes 

15 

800  °C 

600  °C 

2.5 

2  kg  h_1 

Yes 

- 

Yes 

The  height  of  the  biomass  bed  is  detected  by  means  of  the  three 
thermocouples  in  the  reactor.  When  the  gasifier  operates  in  semi- 
continuous  mode,  the  biomass  loads  are  clearly  detected  by  the 
deepest  and  middle  thermocouples.  In  fact,  the  temperatures  reg¬ 
istered  by  those  sensors  inside  the  biomass  bed  are  around  100  °C, 
lower  than  the  gasification  temperature  that  is  externally  imposed. 
In  the  continuous  mode,  the  biomass  is  loaded  and  discharged 
continuously  and  the  height  of  the  bed  is  detected  by  the  middle 
thermocouple. 

3.3.  Catalytic  filter 

In  order  to  directly  feed  a  SOFC  stack  with  the  syngas  obtained 
from  gasification,  a  high  temperature  filtration  system  is  required. 
This  filter  is  required  to  have  high  cleaning  efficiency  at  the  opera¬ 
tion  temperature  (700-800  °C)  and  possibility  to  be  regenerated  or 
easily  disposed;  moreover,  considering  that  gasification  is  carried 
out  at  atmospheric  pressure,  the  filter  must  have  a  low  head  loss. 

The  hot  gas  cleaning  unit  consists  of  a  cylindrical  vessel  of 
60.33  cm  of  diameter  and  40  cm  of  length  and  it  is  placed  down¬ 
stream  the  gasifier.  It  is  insulated  and  heated  up  to  81 5  °C  by  means 
of  an  electric  oven.  The  aim  of  the  catalyst  is  to  reduce  both  the 
tar  concentration,  which  decreases  the  fuel  cell  efficiency,  and  the 
hydrogen  sulphide  that  is  poisonous  to  the  cells.  A  two  stage  fil¬ 
ter  was  tested,  consisting  of  a  bed  of  calcinated  dolomite  used 
for  tar  cracking  followed  by  a  catalyst  made  of  manganese  oxide 
supported  on  zirconium  silicate  for  hydrogen  sulphide  adsorption. 

The  efficiency  of  dolomite  for  tar  cracking  at  high  temperature 
(700-900  °C)  has  been  investigated  and  proved  by  several  authors 
[  1 5-1 7].  On  the  contrary,  for  the  hydrogen  sulphide  abetment,  only 
catalysts  guarantee  up  to  a  temperature  of  450  °C  are  now  com¬ 
mercially  available.  These  catalysts  are  based  on  ZnO  which,  in 
a  reducing  environment,  turns  in  metallic  Zn  that  evaporates  at 


600  °C.  For  process  temperatures  between  600  and  1000°C  sev¬ 
eral  catalysts  are  under  investigation,  but  nothing  is  still  available 
for  industrial  and  commercial  applications.  Manganese  oxide  is 
a  promising  material  for  high  temperature  desulfurization  for  its 
good  performances  and  low  cost  [18]. 

One  of  the  main  drawback  in  hot  gas  cleaning  consists  in  the 
catalysts  long  regeneration  time  due  to  the  low  oxygen  diffusion 
rate  into  the  bulk  of  the  sulphide  catalyst.  Some  authors  proposed  to 
work  with  high  space  velocities  and  short  adsorption  time  in  order 
to  use  only  the  surface  of  the  catalyst,  letting  the  bulk  unreacted 
[19];  this  allows  to  obtain  regeneration  times  similar  to  operation 
times  but,  on  the  other  hand,  leads  to  a  inefficient  utilization  of 
the  catalyst,  being  the  bulk  unused.  A  fine  dispersion  of  sorbent 
material  on  a  porous  support  may  increase  the  catalyst  adsorption 
capacity  with  respect  to  the  sorbent  mass  and,  at  the  same  time, 
make  possible  a  fast  regeneration  process. 

The  catalyst  preparation  follows  the  step  described  below: 


1.  Semi  calcinated  dolomite  Semidol  type  K9/LB  (Franz  Mandt 
Gmbh)  with  a  mean  grain  size  of  2  mm  has  been  calcinated  at 
1000°C  for  4h. 

2.  Zirconium  silicate  supports  have  been  produced  by  dipping 
polyester  open  cell  sponge  cylinders  (d  =  6mm,  L  =  20mm,  pore 
density  =  60  ppi,  Eurotec)  in  a  zirconium  silicate  slurry.  The  slurry 
was  obtained  by  mixing  zirconium  silicate  powder  with  Darvan 
C-N  (RT  Vanderbilt)  and  distilled  water  with  a  weight  ratio  of 
40/1/1 0,  in  Turbula  for  6  h.  After  mixing,  1  g  of  Duramax  B-l 000 
binder  (Rohm  and  Haas)  was  added  for  15.3  g  of  suspension. 

3.  The  supports  have  then  been  dried  in  oven  at  80  °C  overnight 
and  then  calcinated  according  to  the  following  programme:  RT 
to  300  °C  in  5  h,  300  °C  to  1350  °C  in  5  h  50  min,  dwell  at  1350  °C 
for  3  h,  free  cooling. 


Table  5 

Average  gas  composition  of  the  tests  performed. 


Test  no. 

Char  (g  kg-1) 

C2H2  (%  vol) 

C2H4  (%  vol) 

C2H6  (%  vol) 

CH4  (%  vol) 

CO  (%  vol) 

C02  (%  vol) 

H2  (%  vol) 

H2S  (ppm) 

LHV  (MJNm-3) 

1 

n.a. 

0 

0.00 

0.00 

1.3 

4.2 

29.2 

63.4 

- 

7.85 

2 

175 

0 

0.00 

0.00 

1.5 

7.0 

28.2 

63.1 

- 

8.22 

3 

150 

0 

0.00 

0.00 

1.4 

3.7 

30.0 

64.9 

- 

7.96 

4 

n.a. 

0 

0.01 

0.01 

1.9 

6.6 

27.2 

64.3 

- 

8.46 

5 

177 

0 

0.03 

0.03 

2.9 

3.6 

30.4 

63.0 

- 

8.35 

6 

180 

0 

0.03 

0.03 

2.9 

3.6 

30.4 

63.0 

- 

8.32 

7 

8 

9 

200 

0 

0.03 

0.03 

3.5 

4.0 

28.4 

64.0 

- 

8.71 

172 

0 

0.91 

0.40 

6.46 

11.8 

29.1 

51.2 

_ 

10.13 

10 

182 

0 

0.85 

0.37 

6.24 

12.3 

28.15 

52.1 

- 

10.15 

11 

177 

- 

- 

- 

6.37 

9.8 

25.7 

53.72 

84 

9.3 

12 

n.a. 

0.1 

- 

0.6 

6.2 

13 

25.7 

52.0 

15 

9.9 

13 

180 

0.1 

- 

0.44 

7.47 

9.9 

29.2 

52.8 

25 

9.9 
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4.  20  g  of  supports  have  been  introduced  in  a  jar  containing 
lOOmL  of  a  0.5  M  solution  of  Manganese  nitrate  tetrahydrate 
(Sigma-Aldrich)  and  gently  stirred  for  1  h;  after  drying  in  oven 
at  80  °C  for  6  h,  the  supports  were  finally  calcinated  at  850  °C  for 
2  h. 

Before  adding  the  catalyst  to  the  gasifier  filter  unit,  it  has  been 
tested  in  a  alumina  tubular  reactor  heated  in  a  furnace  (Fig.  6).  The 
reactor  can  be  fed  with  different  gases.  Two  tests  have  been  carried 
out:  in  the  first  0.5  g  of  fresh  catalyst  kept  at  600  °C  has  been  fed 
with  1000  ppm  of  H2S  in  H2;  in  a  subsequent  test  the  same  amount 
of  fresh  catalyst  has  been  tested  with  a  simulated  syngas  composed 
by  234  ppm  of  H2S,  3%  of  CH4,  20%  CO,  12%  C02,  20%  H2  and  45%  of 
N2.  The  hydrogen  sulphide  concentration  downstream  the  catalyst 
is  measured  by  means  of  a  portable  Hp  Agilent  3000  Micro  GC. 

The  catalyst  has  shown  good  performances  when  fed  with 
hydrogen  containing  1000  ppm  of  H2S,  being  the  concentration  of 
H2S  at  the  outlet  lower  than  10  ppm  for  40  min.  The  inlet  flow  rate 
varied  from  40  to  200NmLmin_1  so  that  the  amount  of  sulphur 
adsorbed  is  5.8  mgs  g^  (Fig.  7a). 

In  the  experiment  with  the  simulated  syngas  containing  an  inlet 
concentration  of  234  ppm  of  H2S,  the  catalyst  has  adsorbed  H2S  for 
almost  2  h.  Considering  a  gas  flow  of  1 07.4  NmL  min-1 ,  the  amount 
of  sulphur  adsorbed  is  8.6  mgs  g^  (Fig.  7b). 

The  H2S  concentration  in  the  syngas  produced  by  the  steam  gasi¬ 
fier  is  around  85  ppm  (test  no.  11).  The  gasification  test  nos.  12  and 
13  have  been  run  with  the  presence  of  the  catalyst.  The  catalytic 
filter  has  been  filled  with  a  mixture  of  calcined  dolomite  (300  g) 
and  300  g  of  silica  support  impregnated  with  manganese  oxides. 

The  results  are  reported  in  Table  5.  The  catalyst  has  shown 
positive  performances,  however  the  gas  cleaning  session  could  be 
improved  to  guarantee  higher  syngas  quality. 

3.4.  SOFC  performance 

Finally,  the  gasifier  has  been  coupled  with  a  SOFC  stack  and  some 
preliminary  tests  have  been  performed.  Even  if  the  gasifier  has  been 
designed  to  be  coupled  with  a  5  kW  fuel  cells  stack,  the  preliminary 
tests  have  been  run  with  smaller  stacks  of  250  W  (test  no.  14)  and 
330  W  (test  no.  1 5)  to  verify  the  feasibility  of  the  coupling  between 
the  gasifier  and  the  SOFCs  stack. 


Fig.  7.  H2S  abetment  considering  a  mixture  of  H2  and  H2S  concentration  of 
1000  ppm  in  (a)  and  a  simulated  syngas  with  234  ppm  of  H2S  in  (b). 

The  stacks  are  designed  to  have  an  electric  efficiency  of  45%  (on 
the  fuel  UTV)  at  70%  of  the  fuel  utilization.  These  fuel  cell  stacks 
achieves  a  specific  power  density  of  0.3  W  cm-2  of  active  area  at 
a  cell  potential  of  0.7  V  and  a  current  density  of  0.428  A  cm-2.  The 
active  area  of  each  cell  is  50  cm2.  Considering  the  characteristics 
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Fig.  8.  Voltage,  current  and  power  density  of  the  test  no.  14  (power  stack:  250  W). 


of  the  cells,  3  clusters  of  6  cells  each  are  needed  to  reach  250  W, 
and  4  clusters  for  330  W.  Considering  the  enthalpy  of  the  reactions 
(14)-(16)  and  the  moles  of  H2,  CO  and  CH4  in  the  inlet  syngas,  the 
whole  reaction  enthalpy  is  -11,755  kj  per  hour,  corresponding  to 
an  input  power  of  3.26  kW. 

H2  +  l/202^  H20  (14) 

C0  +  1/202^C02  (15) 

CH4  +  202  ->  2H20  +  C02  (16) 

The  number  of  oxygen  moles  required  for  100%  fuel  utilization 
in  the  fuel  cell  is  calculated  with  (17). 

Mol  of  02  =  0.5  (mol  of  H2)  +  0.5  (mol  of  CO)  +  2  (mol  of  CH4)(17) 
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Fig.  9.  Power  produced  and  inlet  and  outlet  fuel  temperature  of  test  no.  15  (power 
stack:  330  W). 


decreased  again  and  the  test  has  been  stopped.  The  cell  efficiency 
has  been  assessed  at  18%  during  the  first  phase  (60  W  produced) 
and  at  40%  during  the  second  phase  (130-140  W  produced). 

It  is  worth  to  highlight  that  the  tested  fuel  cells  are  themselves 
prototypes  of  a  local  industry,  and  some  problems  related  to  the 
non-correct  operation  of  the  SOFCs  can  occur  during  the  gasifica¬ 
tion  tests. 

4.  Results  and  discussion 

4A.  Energy  and  carbon  balance 


Considering  the  average  syngas  composition  and  production 
measured  during  the  experimental  activity,  the  yields  in  mol  per 
hour  are:  H2  =  31.6 mol h-1;  CO  =  7.42 mol h-1;  CH4  =  3.9molh-1 ; 
C02  =  1 7.3  mol  h- 1 ,  for  a  feeding  rate  of  2  kg  h- 1 .  By  means  of  ( 1 7) 
the  moles  of  oxygen  needed  for  the  full  fuel  utilization  have  been 
computed  as  27.34  mol  h_1.  Thus,  the  flow  rate  is  0.66  Nm3  h-1  of 
oxygen,  or  3.2  Nm3  h-1  of  air. 

Since  the  fuel  cell  stacks  tested  are  undersized  with  respect  to 
the  available  syngas,  only  15-20%  of  the  produced  syngas  is  fed 
in  the  stack.  All  the  syngas  generated  at  the  tested  feeding  rate 
(2kgh-1)  can  supply  a  2.5  kW  SOFC  stack,  considering  a  fuel  uti¬ 
lization  of  70%. 

In  Fig.  8  the  voltage,  current  and  power  of  the  fuel  cell  stack 
measured  during  the  test  14  are  reported.  The  stack  has  worked  for 
approximately  80  min  before  stopping  due  to  the  poisoning  of  the 
cells.  There  has  been  a  peak  of  power  production  of  65  W,  after  a 
period  of  increasing  and  decreasing  of  the  power.  Before  the  switch¬ 
ing  off,  an  almost  stable  period  of  operation  has  been  registered 
for  20  min  with  a  power  production  of  30-35  W.  Since  the  nom¬ 
inal  power  is  250  W  (using  pure  hydrogen  as  fuel),  the  stack  has 
worked  at  18%  of  its  nominal  power  (considering  the  performance 
of  the  more  reliable  stable  period).  During  the  initial  peak  of  pro¬ 
duction,  the  stack  has  worked  at  24%  of  its  nominal  power,  even  if 
for  a  short  time. 

In  the  last  test  (no.  1 5)  a  stack  of  24  cells  has  been  coupled  with 
the  gasifier  (4  clusters)  for  a  nominal  power  of  330  W. 

Fig.  9  shows  the  power  generated  by  the  SOFC  stack  and  the  fuel 
input  and  output  temperatures.  In  the  first  phase,  60  W  have  been 
generated  for  approximately  10  min,  then  the  stack  performance 
has  rapidly  decreased.  The  feeding  of  the  syngas  has  been  stopped 
for  1 0  min  to  clean  the  SOFCs  with  the  so  called  “forming  gas”  (nitro¬ 
gen  with  5%  of  hydrogen).  In  the  second  part  of  this  test,  130W 
have  been  produced  for  15  min,  then  the  SOFCs  performance  has 


The  energy  and  carbon  balance  have  been  assessed  on  the  basis 
of  all  the  performed  tests. 

The  carbon  balance  has  been  assessed  using  the  biomass  mass 
flow  rate  and  elemental  composition  as  inputs  and  the  products 
(syngas  and  char)  composition  and  specific  production  per  kg  of 
biomass  as  outputs.  From  the  char  elemental  composition,  reported 
in  Table  6,  it  has  also  been  possible  to  calculate  the  char  heating 
value.  Moreover,  in  order  to  compute  the  carbon  balance  a  value  of 
60  mg  Nm-3  for  tar  production  has  been  considered  [20]. 

The  energy  balance  has  been  also  estimated  considering  the 
input  and  output  energy  fluxes  measured  by  means  of  the 
biomass,  syngas  and  char  characteristics  and  thanks  to  some  energy 
meters  installed  on  the  experimental  plant.  Those  devices  mea¬ 
sure  the  electric  consumption  of  both  the  steam  generator  and  the 
reactor. 

The  carbon  balance  has  been  closed,  for  all  the  tests,  with 
an  error  of  ±8%.  The  system  overall  energy  efficiency  has  been 
assessed  around  40%.  This  value  is  in  agreement  with  the  effi¬ 
ciency  calculated  by  means  of  the  thermodynamic  model  described 
above  (Table  1).  The  energy  required  by  the  electric  ovens  and 
by  the  steam  generator  has  been  compared  with  the  theo¬ 
retical  calculations  obtained  by  means  of  the  thermodynamic 
model.  The  comparison  between  the  predicted  and  measured 


Table  6 

Char  properties. 

Char  ultimate  analysis  (%  mass  a.r.) 

C:  91.35  H:1.75  0:6.73  N:0.17  S<0.01 

Moisture  content  (%  mass  a.r.)  1 .8% 

Ash  (%  mass  a.r.)  1.5% 

Lower  heating  Value  (MJ  kg-1 )  28.7 
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Fig.  10.  Comparison  between  calculated  (lines)  and  measured  (points)  energy  data 
of  (a)  the  steam  generator  at  different  steam  temperature  and  feeding  rate  (b)  the 
reactor  at  different  SC  and  steam  temperature  of  600  °C. 

data  is  reported  in  Fig.  10  and  a  good  agreement  can  be  clearly 
observed. 

5.  Model  calibration 

The  syngas  composition  and  LHV  predicted  by  the  equilibrium 
model  for  the  steam  gasification  process  (SC  =  2,  SC  =  3)  versus  the 
average  gas  composition  measured  is  reported  in  Figs.  11  and  12. 
For  the  first  group  of  tests  the  model  shows  a  good  agreement  for 
the  H2  prediction,  an  overestimation  of  the  CO  and  underestimation 
of  C02.  Anyway,  the  lower  heating  values  predicted  by  the  model 
are  in  good  agreement  with  the  values  calculated  considering  the 
average  syngas  composition  for  the  first  series  of  test.  For  the  sec¬ 
ond  group,  the  higher  methane  concentration  (not  predicted  by  the 
model)  increases  the  syngas  heating  value. 

As  shown  by  previous  literature  studies,  a  correct  methane  esti¬ 
mation  by  means  of  a  thermodynamic  model  is  a  difficult  task,  since 
it  is  not  an  equilibrium  compound  [21  -23  ].  The  high  concentrations 
of  methane  -  measured  in  a  real  reactor  -  can  result  from  incom¬ 
plete  conversion  of  pyrolysis  products  (tars)  and  thus  cannot  be 
explained  on  a  purely  thermodynamic  basis. 

In  the  first  seven  tests  the  modeling  results  are  in  good  agree¬ 
ment,  since  the  methane  production  is  very  low  due  to  the  reactor 
configuration  that  allows  a  high  residence  time.  In  the  second  series 
of  tests,  the  methane  production  is  significantly  higher  with  values 
around  7%. 

Flowever,  the  gas  heating  value  predicted  by  the  model  is  in 
good  agreement  with  the  experimental  values. 

As  a  whole,  many  authors  have  proposed  different  approaches 
in  order  to  achieve  a  better  agreement  between  equilibrium  model 
predictions  and  experimental  data.  A  possible  approach  is  the 
application  of  empirical  parameters  in  order  to  modify  the  carbon 


Fig.  11.  Syngas  composition  foreseen  by  the  equilibrium  model  for  SC  =  2  (a)  and 
SC  =  3  (b)  versus  experimental  average  gas  composition  (points). 


conversion  or  to  correct  directly  the  methane  fraction  in  the  syngas 
[21,22,24]. 

The  high  percentage  of  the  residual  fraction  formed  during  the 
gasification  process  contains  approximately  40%  of  the  initial  moles 
of  carbon.  These  moles,  which  remain  in  solid  phase,  do  not  partic¬ 
ipate  in  carbon  monoxide,  methane  and  carbon  dioxide  formation. 
The  model  has  been  calibrated  according  to  the  experimental  data 
to  take  into  account  the  residual  solid  phase  and  the  methane  pro¬ 
duction.  The  carbon  conversion  efficiency  (i.e.,  the  ratio  between 


Temperature /°C 

— SC=2  —  SC=2.5  -*-SC=3  ♦  SC=3  •  SC=2.5  a  SC=2  ■  SC=3 


Fig.  12.  Syngas  composition  foreseen  by  the  equilibrium  model  for  SC  =  2  (a)  and 
SC  =  3  (b)  versus  experimental  average  gas  composition  (points). 
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H2  CO  C02  CH4 

■  exp  E  original  model  DMOD1  SMOD2 


Fig.  13.  Comparison  between  the  experimental  data  (exp)  and  the  original  and  mod¬ 
ified  model  for  test  no.  9.  MODI  is  the  model  modified  with  the  r]c  parameter,  and 
the  MOD2  includes  also  the  ni,  ri2  parameters. 


the  number  of  carbon  moles  converted  in  gas  and  the  total  carbon 
moles  fed)  has  been  calculated  with  (18). 


(18) 


Then  the  initial  input  vector  of  the  model,  which  contains  the 
number  of  moles  of  C,  H,  O,  N,  S  per  kg  of  biomass  fed,  has  been 
modified  according  to  the  value  (19). 

Kput  =  ["c,  nH,  no]  =►  Nlnput  =  [rjcric,  nH>  n0]  (19) 


— H2  -o-CO  —C02  — CH4  — H2S 

Fig.  14.  (a)  Reactor  temperature  simulated  by  the  finite  element  model;  (b)  syngas 
composition  simulated  by  the  modified  equilibrium  model  for  z  =  456  mm. 


To  take  into  account  also  the  methane  formation,  the  experi¬ 
mental  gas  composition  has  been  used  to  evaluate  the  moles  of 
carbon  (ni )  and  hydrogen  (n2)  converted  into  methane  during  the 
process.  The  initial  composition  is  corrected  considering  the  r)c  and 
subtracting  the  moles  of  carbon  and  hydrogen  arising  from  the 
previous  calculation  (20). 

NinPut  =  focnc,  nH,  n0]  =>  N*nput  =  [/jcnc  -nunH-  n2,  n0]  (20) 

In  Fig.  1 3  the  experimental  gas  composition  of  test  no.  9  has  been 
compared  with  the  original  and  the  modified  models  outputs. 

The  agreement  between  experimental  and  modeling  data,  con¬ 
sidering  only  the  carbon  conversion  parameter  (output  of  MODI), 
has  been  remarkably  improved.  However,  in  order  to  achieve  a  bet¬ 
ter  agreement,  is  worth  to  also  calibrate  the  methane  concentration 
in  the  syngas,  as  clearly  shown  by  the  output  of  MOD2.  The  range 
of  temperature  and  SC  investigated  in  the  present  experimental 
activity  is  not  wide  enough  to  find  out  experimental  correlations 
capable  to  estimate  the  methane  and  carbon  conversion.  However, 
as  found  by  other  authors  through  extensive  experimental  cam¬ 
paigns,  it  is  possible  to  find  out  experimental  correlation  between 
the  gasification  parameters  and  carbon  conversion,  i.e.,  between 
the  carbon  conversion  and  the  ER  value  [18].  From  one  side  the 
use  of  empirical  correlations  limits  the  predictive  capability  of  a 
thermodynamic  model  to  specific  process  conditions  (i.e.,  reac¬ 
tor  type  and  design).  However,  from  the  other  side,  the  resulting 
‘quasi-equilibrium’  model  becomes  an  useful  tool  to  assess  the  gas 
composition  and  yield  with  satisfying  accuracy  for  configurations 
similar  to  the  test  case. 

To  test  the  reliability  of  the  2D  FEM  model  described  in  Sec¬ 
tion  2.3,  the  temperature  data  measured  by  the  thermocouples 
placed  inside  the  reactor,  at  different  radial  positions  and  heights, 
have  been  compared  with  the  temperature  profiles  simulated  by 


the  model.  The  experimental  temperatures  range  between  700  and 
800  °C,  while  the  computed  ones  are  between  500  and  600  °C. 

The  finite  element  model  takes  into  account  only  the  heat  trans¬ 
fer  through  conduction.  The  deviations  between  the  experimental 
data  and  the  model  outputs,  have  confirmed  the  influence  of  the 
heat  exchanges  through  radiation  and  convection.  In  order  to  cal¬ 
ibrate  the  model  response,  a  single  heat  transfer  coefficient  (i.e., 
thermal  transmittance)  has  been  estimated  using  the  experimen¬ 
tal  data.  The  estimated  thermal  transmittance  has  been  defined 
as  an  “apparent  thermal  conductivity”  including  the  contribu¬ 
tion  of  the  convection  and  radiation.  For  this  purpose,  the  Fourier 
equation  for  heat  conduction  (21)  in  one-dimensional  case,  using 
non-dimensional  parameter  (22)  has  been  considered. 


d2o  n 

dy2  +2y  dy 


=  0 


(21) 


T-Tj  _  _  X 

To-Tj'  Y~  2Var 


(22) 


A  time  dependent  boundary  condition  has  been  applied,  in 
agreement  with  the  heating  rate  of  the  oven.  The  solution  of  heat 
equation  in  a  semi-infinite  solid  body  has  been  assumed  for  the 
purpose  of  the  present  analysis.  The  solution  is  represented  by  the 
error  function  (erf). 

6)=^3^  =  1-^  jV  e-z2dz  =  \-erf{y)  (23) 

The  solution  of  the  integral  of  the  erf  function  is  not  known  in 
an  explicit  form.  The  values  of  the  error  function  are  tabled  [25]. 
Thanks  to  the  temperature  measurements,  it  has  been  possible  to 
obtain  the  heating  time  for  the  biomass  starting  from  the  initial  7) 
to  the  final  temperature  Tf.  Thus,  the  <9  and  the  y  values  have  been 
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calculated.  From  the  definition  of  y  it  is  possible  to  deduce  A,  the 
apparent  thermal  conductivity  of  the  solid  body  (24). 


The  apparent  thermal  conductivity  adopted  in  the  finite  element 
model  is  an  average  of  the  A  values  calculated  by  means  of  the 
experimental  data  of  the  experimental  tests.  The  estimated  thermal 
conductivity  is  1 .23  W  rrr 1  K-1 . 

The  2D  model  -  corrected  with  the  apparent  thermal  conductiv¬ 
ity  -  has  been  coupled  with  the  quasi-equilibrium  model,  to  take 
into  account  the  biomass  conversion  in  char  and  methane.  The  tem¬ 
perature  profiles  are  shown  in  Fig.  14.  One  of  the  thermocouples  is 
located  at  235  mm  from  the  top  of  the  reactor  and  equidistant  from 
the  reactor  wall  and  the  central  coaxial  tube.  The  data  measured 
by  the  sensor  show  a  good  agreement  with  the  simulated  temper¬ 
ature.  The  predicted  temperature  profiles  inside  the  reactor  can  be 
reasonably  considered  close  to  the  real  one. 

6.  Conclusions 

At  present,  the  biomass  thermochemical  conversion  processes 
and  their  applications  for  heat  and  power  generation  are  object  of 
several  research  studies,  to  find  out  efficient  systems  to  exploit  the 
energy  content  of  biomass. 

The  present  work  aims  to  improve  the  knowledge  on  the  poten¬ 
tial  of  the  steam  gasification  process  for  power  generation  in  small 
scale  applications.  The  syngas  produced  via  steam  gasification  pro¬ 
cess  seems  to  be  a  suitable  fuel  both  for  internal  combustion 
engines  and  for  solid  oxide  fuel  cells  due  to  its  high  hydrogen  con¬ 
tent. 

Engineering  models  are  usually  applied  to  design  and  evalu¬ 
ate  the  performances  of  gasification  systems.  A  thermodynamic 
equilibrium  model  has  been  used  to  estimate  the  syngas  compo¬ 
sition.  The  positive  aspect  of  the  thermodynamic  approach  lays 
in  its  applicability  to  several  systems  without  a  deep  knowledge 
of  the  reaction  mechanisms.  Thus,  it  can  be  successfully  applied 
to  determine  the  maximum  theoretical  performance  of  a  biomass 
conversion  process. 

A  continuous  small  scale  steam  gasification  system  has  been  set 
up  and  tested.  The  syngas  produced  has  a  hydrogen  content  rang¬ 
ing  between  50  and  60%  and  a  UTV  of  8MJNm-3.  Considering  its 
composition  and  energy  content,  the  obtained  syngas  is  a  suitable 
fuel  for  fuel  cells.  However,  the  gas  cleaning  is  still  one  of  the  main 
critical  issues.  In  particular  the  tar  and  the  H2S  in  the  gas  can  rapidly 
decrease  the  life  of  the  fuel  cells.  The  tested  catalyst  shows  a  satis¬ 
fying  efficiency  in  H2S  and  tars  abatement  even  if  for  a  short  period 
of  operation,  after  which  the  stack  performances  decrease.  This  fact 
may  be  attributed  to  filter  undersize. 

A  satisfactory  agreement  has  been  found  between  the  experi¬ 
mental  data  and  the  modeling  simulations  calibrated  by  means  of 
the  data  collected  during  the  experimental  activity.  A  good  agree¬ 
ment  has  been  also  observed  between  the  calculated  and  measured 
energy  fluxes  required  by  the  system  in  order  to  sustain  the  con¬ 
version  process. 


The  thermodynamic  approach  has  been  confirmed  to  be  a  simple 
engineering  tool  useful  to  assess  with  good  reliability  the  theoret¬ 
ical  performance  of  a  gasification  system,  knowing  only  the  main 
gasification  parameters. 
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